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Abstract 

We investigate the percentage of delivering capacities that are actually consumed in a typical 
traffic dynamics where the capacities are uniformly assigned over a scale-free network. Theoretical 
analysis, as well as simulations, reveal that there are a large number of idle nodes under both free 
and weak congested state of the network. It is worth noting that there is a critical value of effective 
betweenness to classify nodes in the congested state, below which the node has a constant queue 
size but above which the queue size increases with time. We also show that the consumption ratio 
of delivering capacities can be boosted to nearly 100% by adopting a proper distribution of the 
capacities, which at the same time enhances the network efficiency to the maximum for the current 
routing strategy. 
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I. INTRODUCTION 



The interplay between the dynamics and the topology of network can be used to model 
the evolutions and structures of a wide variety of complex systems in nature and human 
society. A branch of dynamics research that focuses on the data traffic in the network has 
drawn a lot of interests. The studies in this field are mainly about capturing the attributes 
of the data transportation and providing hints on how to make the physical communica- 
tion network more efficient. Previous structure analysis based on real data spanning several 
disciplines revealed that many networks, including the Internet, exhibit scale-free character- 
istics 111. The degree distribution of these networks follows a power law. Having explored 

nn 

the topological properties |2J, |3J , researchers habitually choose the scale-free network as the 
underlying structure when studying the traffic dynamics in order to observe phenomena that 
are close to those in the real world 4 , 0, Q . 

The traffic models proposed by recent works generally consist of four components: the 



underlying network topology, the information packet generation, the routing strategy and 

i, ^ | . They all have counterparts 



the packet delivering capacity of individual nodes [l, 
in the reality. Take the Internet as an example. The four components correspond to the 
physical links between computers, the average data flow, the algorithm implemented on 
routers and the number of parallel packet processors per site, respectively. They also have 



12| and organization design 
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In 



other interpretations in the field of path navigation 
these models, a continuous phase transition of network from free state to congested state 
is observed when the packet generation rate exceeds a critical value [io|. The core purpose 
of researches on these models is the optimization that enhances the network capacity which 
is measured by the critical value of packet generation rate. Some of the previous efforts 
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while the others are 
However, in most of the 



concentrate on finding a most suitable network topology \A 
mainly about deciding an optimal routing strategy 
models, constant delivering capacity is assigned to nodes. It is not economical as the loads 
on different nodes vary. In fact, the delivering capacities are processor resources that they 
should also be taken into consideration during optimization attempts. 

With the hope of discovering potential possibilities of improvements on the communi- 
cation network, we analyze the percentage of delivering capacities that are truly active in 
our traffic model where the capacities are uniformly distributed over a scale free network. 
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It is concerned that congestion is sometimes inevitable in real lives, for example, everyone 
has the experience of traffic jam in rush hours, so the performance of delivering processors 
under congested state of network is paid special attentions. We find that some of the nodes 
in a weak jammed network always have constant number of packets waiting to be delivered 
which is the same situation as in the free state of the network, whereas the others have 
increasing number of packets. These two types can be distinguished by referring to the 
effective betweenness which is an attribute of node defined in Ref. ^]. Our theoretical esti- 
mates and simulations reveal that delivering capabilities of nodes are not well utilized even 
under the congested state of network, though optimal routing strategy based on local degree 
information discussed in Ref. is adopted. To alleviate the waste of processor resources, 
we apply a non- uniform distribution of the delivering capacities which adjusts the network 
to the best performance under the given conditions. 

This paper is organized as follows. The traffic model is described in Sec. II. In Sec. Ill, 
theoretical analysis and simulations of the model are provided in both free and congested 
state. The conclusion is given in Sec. IV. 

II. TRAFFIC MODEL 

We use the famous model proposed by Barabasi and Albert to build the underlying 
scale-free network. BA model, which features growth and preferential attachment, uncovers 
for the first time the mechanism controlling the emergence of power-law degree distribution 
observed in real networks. We set the model parameters mo = m = 5 and network size 
N = 1000. In the traffic dynamics, network nodes are thought to be sources that produce 
information packets, routers that deliver packets and first-in-first-out queues with unlimited 
size that store to-be-processed packets. At each time step, for a node %{% = 1, . . . , N), the 
following procedures are done. A packet is generated with probability p«. The new packet 
whose destination is randomly selected among the other N — 1 nodes is placed at the back 
of queue i. Meanwhile, packets are picked from the front of queue i. If the picked packet has 
destination node i, it is removed; otherwise, it is delivered to node j, one of the neighbors 
of node i, with preferential weight 
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where sum runs over the neighbors of node i and ki is the degree. In Eq. (JI}, we have 
used the optimal routing strategy based on the knowledge of only neighbor's degree [ll| . At 
most Cj packets can be processed by node i if there are adequate ones in queue i. is the 
delivering capacity of node i. We set pi = p, Ci = C = 10 following the parameters used in 



Ref. in order to meet the conditions for the optimal routing strategy. These procedures 
are carried out for different i in parallel. What we concern is the ratio between the number 
of packets actually processed by the network during a unit time and the total capacities 
assigned, which is defined as the effective delivering capacity. 

III. EFFECTIVE DELIVERING CAPACITY 

De pen ding on p, there is a continuous phase transition from free state to congested 
state [l2j. When p is small, the packets flow freely in the network as the nodes always 
have available processing abilities to send them to the next positions. The average travel 
time t of the packets keeps constant and the total number N p (t) of packets floating in the 
network at time step t fluctuates slightly around pNr. However, when p is large, packets 
beyond the delivering capacities accumulate continuously in the queues, causing r to diverge. 
Since the number of packets the network can manage and the generation rate p are both 
constants, the amount of accumulation at every time step is also constant. As a result, N p (t) 
increases linearly with time. To watch the transition accurately, we use the order parameter 



£ introduced in Ref. 
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where (. . .) indicates an average over time windows of width At. From the properties of 
N p (t), we know that £ is zero in free state but non-zero in congested state. p c , the critical 
value of p, which is also the maximum generation rate keeping the system in free state, 
measures the network capacity. 

To characterize the situation in free state, we use the method introduced in Ref. [8j] with 
our modifications. Let us focus on a packet at node i whose destination is node k. The 
probability for this packet to go to node j in one time step is denoted as p^. The precise 
form of "p\a depends on the routing strategy. For the scheme used in the model, 



where Aij is the element of the adjacent matrix. The strategy is Markovian that each 
packet is delivered independently. As the network is in free state, packets pass through 
nodes without any wait. The probability for the above packet to go to node j in n time 
steps is given by 

4( n )= £ PiiAh-'-PLir ( 4 ) 

h,fa,...,l„-l 

Treating and -P^(n) as elements of matrices p k and P k (n) respectively, we have 

P k (n) = (p fc ) n . (5) 

Let us switch the focus to the centrality of node. In the traffic model, packets coming from 
neighbors at a specific time step are a fraction of packets generated before. Supposing that 
one packet with target node k is generated at node % each time step, we calculate b k j(t), the 
average number of such packets moving to node j at time step t, 

*>*(*) = S> fc ) B - (6) 

n=l 

As N p (t) is stable, when t — > oo, 

oo 

b fe = ^(p fc ) n = (i- P fe )-y\ (7) 

n=l 

Summing over all the possible sources and targets of packets yields the effective betweenness 
of node j, Bj, 

Bj = Y, b ir (8) 

i,k 

Note that Bj depends on both the routing strategy and the network topology. As the routing 
algorithm can vary, highly connected nodes do not always have large effective betweenness. 
The load Lj, which is the average number of packets to be delivered by node j every time 
step, is written as, 

where the first term on the right side corresponds to packets generated by j itself while the 
second term corresponds to packets coming from neighbors. Any node j with Lj > C will 
cause the network to be congested. Thus we estimate p c by 

o ~ ° - C(N - 1] (10) 

N-1 ^ 
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where B* is the maximum effective betweenness. In our model, we had p c = 0.0045. When 
p < p c , the number of packets processed at node j per time step is exactly its load Lj. 
Therefore, effective delivering capacity rj under free state is 

which is proportional to p. 

For p larger than p c , Eq. (JTTJ) is invalid because the system is not expected to be in free 
state when t — > oo in Eq. (JZj). Back to time step t — 1, the average load in the network is 
equal to the generation rate p. Here we limit our discussion to a weak congested scenario 
where p is not supposed to be extremely large so that the nodes are still able to process all 
the packets in the queues for the first few time steps. During these time steps, we can obtain 
the instantaneous load Lj(t) by replacing Eq. (J7J) with Eq. (jUJ) in the previous calculation for 
free state. As Lj(t) is a monotonically increasing function of time, it is at a later time step t c 
that the largest Lj(t) starts to exceed the delivering capacity C . While the packets beyond 
the capacity are hindered at the corresponding node, they do not flow in the network. 
Comparing with an imaginary case where no constraints on the delivering capacities are 
applied, the actual instantaneous load increases with time more slowly. Though Lj(t) for 
t > t c can not be directly worked out by analytical calculation, we can predict that with 
more nodes suffering from packet overloading, Lj(t) increases even more slowly. After a 
transient time, the instantaneous load over the whole network stops increasing and remains 
stationary. Thereafter, all the nodes in the network maintain their states. Those with load 
higher than delivering capacity have growing queues, while the others keep constant queues. 
Since the effective betweenness Bj determines free state load Lj by Eq. (JHJ), it is plausible to 
suppose that nodes with greater effective betweenness become overloaded more easily. We 
expect a critical effective betweenness B c for current p, which can be used as the criterion 
to classify the two types of node. Once B c is found, effective delivering capacity can be 
considered by separate examinations on different types of node. 

To find B c , we look at the congested state of system when t — > oo. The increasing queues 
are already considerably long. A newly hindered packet at a specific node will wait in the 
queue for so long a time that it seems to be disappeared from the views of the other nodes in 
the network. We assume that a node with an increasing queue is a "target" for information 
transportation where packets disappear. To make it more like a "target", the C packets it 
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lets out at every time step are shifted to the neighbors as equivalent generation rates. By 
this imagination, the original congested network turns out to be a smaller one in free state, 
and method for p < p c can be exploited. We use the following algorithm: 

1. All the nodes in the network are sorted by effective betweenness Bj in descending 
order. The first node in the sequence is added to an initially empty set S. 

2. The new generation rate //• for node j is 

f/ j = P + J E c Pij> ( 12 ) 

ies 

where sum runs over all the elements in S and pij represents the routing strategy, 

3. In order to suppress packets arriving at "targets", row i in the original matrix p k is 
set to for each i G S. 

4. Employing p\ and the modified p fc , the new load Lj is evaluated by 

L j={ ( 14 ) 

[P+JTi Ei,fe p'Mj otherwise 

5. The validity of 

max(L'j) < C (15) 

must be ensured, which corresponds to constant queue size for nodes not in S. If there 
exists load L'- larger than C, the next node in the sequence descending ordered by Bj 
is added to S. 

6. Step |21 to are repeated until Eq. (JT5|l is fulfilled. B c for current p is just the critical 
Bj separating nodes within S from those outside. 

To describe the changes of queues in the simulation, we define Kj for node j, 

= Um (Qj(t + At) - Qj{i)) 

3 t->oo At ^ 1 

where Qj(t) is the queue size of node j at time step t and (. . .) indicates an average over 
time windows of width At. Kj is zero if node j has a constant queue or positive if node j has 




FIG. 1: Simulation results of k versus B. The critical point B c is found by calculation. In the 
inset, we plot Hj = (Qj(2t) — Qj{t))/Qj{2t) when t — > oo for nodes with Bj > B c , which shows 
that the queues increase linearly with time. 

an increasing queue. Taking p = 0.02 as an example, we illustrate Kj versus Bj in Fig. 
where B c is obtained by calculation. It is shown that below B c nodes have constant queues 
but above B c the queues increase linearly with time. Besides, Kj ~ Bj — B c for Bj > B c . 

We are ready to estimate the effective delivering capacity under congested state. For the 
nodes with constant queues, the packets delivered are equal to the load L'j. For the nodes 
with increasing queues, they work at the capacity C. Note that during the calculation, 
packets processed by these nodes are shifted to the neighbors, which indicates that the 
packets pass through nodes twice. The effective delivering capacity for congested state, rf , 
is then given as 

= H,L> + 2CN, 
1 CN • ( ' 

where N s is the number of elements in S. 

As shown in Fig. we calculated the effective delivering capacity for different p ranging 
from 0.001 to 0.03 and compared the results with simulation where packets processed in a 
unit time are directly counted. In the free state, r\ increases with p according to Eq. (fTT| . 
In the congested state, though the increment of p shows a trend to enlarge L'j, there is a 
cutoff imposed by the delivering capacity C. With the overloaded nodes ignored, the sum 
of L'j does not vary too much. 7/ increases only due to the increment of N s which is small 
when p changes hardly. 
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FIG. 2: Effective delivering capacity r\ and order parameter £ for C = 10. p c follows Eq. (|1U|) . 

It is undesirable that up to p = 0.03 ~ 6.7p c the effective delivering capacity is no more 
than 75%, despite the use of optimal local routing algorithm. In other words, while the 
network is not able to handle all the information packets, 25% processors are idle. An 
optimization problem arises naturally: given the fixed number of processors, how to exploit 
them all? The key is to avoid resource wastes on less loaded nodes. We introduce a new 
delivering capacity distribution, 

C[ = ^rCN, (18) 

where the total capacity CN is conserved. With respect to Eq. Q, when p is approaching 
the critical value, the delivering capacity just meets the load, and no waste happens. Above 
the critical value of p, cutoff imposed by delivering capacity does not reduce the load on 
the nodes with constant queues, and the utilization ratio keeps 100%. In particular, the 
network capacity is enhanced to the maximum since all the processors are active. Similar 
to Eq. (fTo"|) . for an arbitrary distribution of delivering capacity, the critical point p c is 

Pc = min{ Cl{N ~ l \ j = 1, . . . , N}, (19) 

with constraint J2j Cj = CN. The optimal value of p c is obtained by neutralizing the 
differences among Cj(N — 1)/Bj under the constraint. It can be accomplished by employing 
Eq. (JTHJ). The corresponding critical point p c opt is 

CN(N - 1) 

Pco P t= 1 „ ■ (20) 
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FIG. 3: Effective delivering capacity r/ and order parameter £ for new distribution of C captured 
by simulation. p c op t follows Eq. (I2U1) . 

Calculation yields p c opt = 0.0074. In the simulation, delivering capacity given by Eq. (fTHJ) 
is converted to integer where errors are involved. Figure El illustrates simulation results of 
the new distribution. It is clear that the effective delivering capacity is very close to 100% 
for p > p c opt and p c opt is consistent with theoretical estimates. 



IV. CONCLUSION 



In this paper, we study the effective delivering capacity for a typical traffic model un- 
der both free and weak congested state. It is shown that with equal number of packet 
processors assigned to nodes, a large percentage of the processors are idle even under the 
congested state. By redistributing the processors according to effective betweenness, all of 
them are activated, which corresponds to maximum network capacity for the routing strat- 
egy adopted. Considering the difficulties to change the topologies of real networks, previous 
works P, 0| on optimization mainly focus on finding the routing algorithm best fits the 
underlying structure. We rise the importance of choosing a proper distribution of the de- 
livering capacity which is an essential supplement to the optimized routing strategies. In 
addition, we successfully separate nodes with constant queues from those with increasing 
queues in a weak congested network by referring to a critical value of effective betweenness, 
which may provide some hints for further studies on the congested network. 
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